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ABSTRACT 

We analyse the phase diagram of the lattice gauged Nambu-Jona Lasinio 
model with the help of a mean field approximation plus numerical simu- 
lations. We find a phase transition line in the coupling parameters space 
separating the chirally broken phase from the symmetric phase, which is in 
good qualitative agreement with results obtained in the quenched-ladder ap- 
proximation. The mean field approximation relates the critical exponents 
along the continuous phase transition line with the mass dependence of the 
chiral condensate in the Coulomb phase of standard noncompact QED. Our 
numerical results for noncompact QED strongly suggest non mean field ex- 
ponents along the critical line. 
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The gauged Nambu- Jona Lasinio (GNJL) model has become increasingly 
interesting in recent time, one of the reasons for this increasing interest being 
the possibility to define a strongly coupled QED with non trivial dynamics 
[0]. In fact if a non gaussian fixed point exists in non compact QED, the 
naive dimensional analysis does not applies. Therefore operators of dimen- 
sion higher than four, which are non renormalizable in perturbation theory, 
could acquire anomalous dimensions and become renormalizable @]. A good 
candidate is the four Fermi interaction, which when added to the standard 
QED lagrangian preserving the continuous chiral symmetry, gives us the 
GNJL model. 

The lattice action for the GNJL model with noncompact gauge fields 
and staggered fermions reads 



^ = ^ J2 ^li^i'^) + + - GY,XnXnXn+,,Xn+^^■ (l) 

where 9^^, is the standard noncompact plaquette action, P the inverse square 
coupling, A{6) the massless Dirac operator for Kogut-Susskind fermions and 
G the four fermion coupling. 

In the chiral limit, m = 0, this action is invariant under the continuous 
transformations 

An ' An c Xn ' An c y^J 

which define a continuous chiral U{1) symmetry group. 

The vacuum expectation value of the chiral condensate is given by the 
following ratio of path integrals over the Grassmann and gauge fields 



/ [dOdxdx] e ^ yJ2n XnXn 
{XX) = ri^nj.-.j^ (3) 



-s J. 
v_ 

I [dOdxdx] e-^ 

The main technical difficulty when computing vacuum averages as (3) in 
the GNJL model comes from the fact that the action (1) is not a bilinear 
of the fermion fields. The standard procedure consists in the introduction of 
an auxiliary vector field which allows to bilinearize the fermion action. The 
prize to pay for that is that we have one more field to include in the nu- 
merical simulations of this model that besides the number of free parameters 

m, G), makes it difficult to analyse this model with reasonable computer 
resources [Bl. 



2 



Alternatively, we can perform a standard mean field approximation which 
also bilinearizes the action (1). Following the mean field technique, we make 
in (1) the following substitution 

Gj2XnXnXn+nXn+^i ^ 2ciG( XX ) XnXn (4) 

n,fi n 

where d is the space-time dimension. The action (1) becomes in this way 
a bilinear in the fermion fields and the path integral over the Grassmann 
variables can be done by means of the Matthews-Salam formula. 

The v.e.v. of the chiral condensate (3) after the substitution of the mean 
field approximation (4) in the action (1) is given by 



(xx) 



/ [de] e-ii:e^.(n) ^^^[^ + _ 8G{xx)mtr^^^^^^^l^^ 



J [de] e-i E W det[A + {m - 8G{xx))I] 
which after simple algebraic operations can be written as 



(5) 



1 1 



\ V Pi ^'j + {m- SG{ XX )r I 

where the sum in (6) runs over all positive eigenvalues of the massless Dirac 
operator and the integration measure in the v.e.v. includes the fermionic 
determinant of standard noncompact QED, evaluated at the effective mass 
rh = m — SG{xx )• the chiral limit m = 0, equation (6) becomes 

(XX) = 16G(XX)^^E,.^64^.^-^^.; (7) 



1. The phase diagram 

Equation (7) is always verified if ( xx ) = 0) ^-^d this is the only solution 
in the symmetric phase. In the broken phase where ( xx ) 7^ 0) the v.e.v. of 
the chiral condensate will be given by the solution of the following equation 

/ 1 V/2 ^ y 

which gives for the critical line, where the chiral condensate vanishes contin- 
uously, the following expression 
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GJP) = 3_ (9) 

V \ ^i=l / 

The existence of this critical hne in the GNJL model was discovered 
some time ago IQ in the continuum formulation using the quenched-ladder 
approximation. 

Let us discuss qualitatively the phase diagram. For P < P^, where is 
the critical coupling at G = 0, the symmetry is always spontaneously broken 
since in this case equation (8) has a non vanishing solution for any G ^ 0. 
On the other side, the symmetry is also spontaneously broken in the G ^ oo 
limit since in this limit equation (8) can be written as 

from which it follows that 

(XX>~^ (11) 

hi (3 = oo the theory can be solved analytically in this approximation. 
The value of the critical four fermion coupling is in this limit Gc = 0.2017. 
For G values smaller than this value, the symmetry is restored. 

In Fig. 1 we present our numerical results for the phase diagram in the 
13, G plane. The critical line has been obtained by computing numerically 
the v.e.v. of the sum of the inverse square eigenvalues (eq. (9)), which is 
proportional to the chiral transverse susceptibility of the standard noncom- 
pact QED in the chiral limit. The numerical simulations where performed 
using the MFA approach 0, which allows to do computations in the chiral 
limit. We refer the interested reader to the extended bibliography on this 
subject 1^ and especially to the ref. |0 where the computation of the chiral 
susceptibility and the determination of the critical coupling in noncompact 
QED is discussed in detail. 



2. The critical exponents 

The phase diagram of Fig. 1 is in good qualitative agreement with the 
corresponding phase diagram obtained in the quenched-ladder approximation 
[^]. Using this analytical approach, a line of critical points with continuously 
varying critical exponents was found in 0, the intersection point of this line 
with the G = axis corresponding to an essential singularity 0. 
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Later on, numerical simulations of noncompact QED disproved the es- 
sential singularity behavior [Q, putting in evidence the limitations of the 
quenched-ladder approximation. Since our approach contains weaker approx- 
imations, we do hope to get more reliable results for the critical exponents. 

In order to extract the critical exponents, we will start from the key 
equation of state (eq. (6)) relating the order parameter with the "external 
magnetic field" m and the gauge and four fermion couplings. Using the 
previous notation we can write equation (6) as 

(xx) = -2mF(/3,m) (12) 

where the right hand side in (12) is just the chiral condensate in full non- 
compact QED evaluated at the gauge coupling value /3 and fermion mass 
rh. Concerning critical exponents the interesting physical region, as follows 
from the phase diagram of Fig. 1, is /3 > (3^ (Coulomb phase of noncompact 
QED). 

Since we are interested in the critical region (m — 0, ( xX ) — ^ 0), we will 
analyze the behavior of F{(3,fh) in the m — limit. In this limit we can 
write 

F(/?, m) = F(/3, 0) + BrrT + ... (13) 

The second term in (13) possibly contains also logarithmic contribu- 
tions and F{f3, 0) is half the massless transverse susceptibility in noncompact 
QED. Therefore we can write 

^(/?,0) = ^{Y.^/\))p^^^, = l/16a(/3) (14) 

where Gc{l3) in (14) stands for a generic point of the critical line in Fig. 
1. Equation (13), after the substitution of m by m — 8G'(xx)) implies the 
following behavior for the chiral condensate in the m ^ limit 

{xx)^m^ (15) 
and therefore the uj and 5 exponents are related by the equation 

5 = uj + l (16) 

A straightforward calculation allows to compute also the magnetic I3m 
and susceptibility 7 exponents, the final result being 
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(17) 
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The hyperscaling relation 7 = — 1) is verified, as follows from (17). 

The determination of the critical exponents of the order parameter in our 
mean field approach reduces therefore to the determination of the u exponent 
which controls the mass dependence of the chiral condensate in the Coulomb 
phase of noncompact QED. In the /3 — 00 limit of noncompact QED, the 
theory is free and the chiral condensate can be analytically computed. The 
well known result in this case {u = 2 plus logarithmic corrections) implies 
mean field exponents for the end point of the phase transition line, with the 
following behavior for ( xx )i3=oo,g=g^ 



In the general case, the chiral condensate in the Coulomb phase of non- 
compact QED {{xx) ncqed) can be parameterized as follows 



The first contribution in (19) is linear in m, as follows from the fact 
that the massless transverse susceptibility is finite in the Coulomb phase 
of noncompact QED. The next contribution can possibly have logarithmic 
corrections, as happens in the (3 ^ 00 limit where it becomes m^logm. In 
order to extract the uj exponent from the numerical simulations, we can use 
the results for the massless chiral transverse susceptibility to fix A{j3) in 
(19) and fit the numerical results with eq. (19). This procedure has the 
inconvenient that higher order contributions in (19) can induce systematic 
errors in the determination of cu. A better strategy is to measure the massless 
nonlinear susceptibility, defined as the third mass derivative of the chiral 
condensate. In this case we get only one contribution in the m — > limit 
which is logarithmically divergent in the free field theory against a power 
divergence, which will appear if u; < 2. 

Of course in a finite lattice, the non linear susceptibility is always finite. 
However simple finite size scaling arguments tell us that the nonlinear sus- 
ceptibility should diverge logarithmically with the lattice size in the free field 
case whereas a power divergence with the lattice size is expected in the case 
uo <2. In Fig. 2 we have plotted our results for the nonlinear susceptibility 
Xni of noncompact QED against the lattice size aX (3 = 0.237, a value which 
is unambiguously in the Coulomb phase of this model |]l|. This is a log- log 




(18) 



{xx)NCQED = A{[3)m + B{l3)rrf^^ + ... 



(19) 
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plot and the four points correspond to lattice sizes 4, 6,8 and 10. As it is 
shown in the figure, the four points are very well fitted by a straight line, 
this implying that lo <2. 

Due to the potentialities of the MFA method, we have computed vacuum 
expectation values of other operators, which can be considered as generaliza- 
tions of a term contributing to the masslcss nonlinear susceptibility. More 
precisely we have defined Xg by the expression 




When g = 4, we get one of the contributions to the standard massless 
nonlinear susceptibility. In the general case we can write this vacuum ex- 
pectation value as an integral over the spectral density of eigenvalues in the 
following way 




and if the density of eigenvalues p(A) behaves like near the origin, Xq will 
diverge when q > p + 1. In such a case and for lattices of finite size, we 
expect for Xq the following behaviour with the lattice size L 

Xg ~ (22) 

where a in (22) is some positive number. 

It is interesting to note that the p exponent which controls the small A 
behavior of the spectral density p(A), can be related to the cu exponent by 
the following equations 

LU = p — l{p < 3) 

u = 2{p> 3) (23) 

These relations allow to extract the u exponent from the finite size be- 
havior of the generalized nonlinear susceptibility Xg- 

In Fig. 3 we have plotted our results for the inverse of the generalized 
nonlinear susceptibility Xg against the inverse lattice size for q values running 
from 2 to 4 and (3 = 0.237. The solid lines in this figure correspond to a fit 
of all the points at any fixed q with the function 
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X-\L) = a, + \L-' (24) 

The results reported in this figure show that, in the infinite volume limit, 
the inverse generalized nonlinear susceptibility vanishes at large q and is 
different from zero at small g, as expected. Fig. 4 is a plot of the extrapolated 
values of Xq^ (thermodynamical limit) against q. The critical value of q 
at which Xq^ vanishes can be estimated from these results. Hence we get 
qc ~ 2.5 at /3 = 0.237, which implies p ~ 1.5 and uj ~ 0.5. Using now 
the relations (16), (17) the following results for the order parameter critical 
exponents can be derived 

5-1.5, /5„~2, 7=1, (25) 

values which are clearly outside the range of the mean field exponents. 

Non mean filed exponents within a mean field approximation might seem 
rather surprising at first sight. There is however no real contraddiction, since 
we have applied the mean field approximation to the fermion field, while 
fiuctuations of the gauge field are fully taken into account in our numerical 
simulations. In the infinite (3 limit, where the gauge field is frozen to the free 
field configuration, we get mean field exponents. However fiuctuations of the 
gauge field at finite (3 seem to play a fundamental role in driving critical 
exponents to non mean field values. 

The picture which emerges from this calculation is that the critical expo- 
nents change continuously along the critical line of Fig. 1 from their mean 
field values (end point of the critical line) to some (non mean field) values at 
the critical point of noncompact QED. The 6 exponent approaches its mean 
field value [5 = 3) from below whereas the magnetic exponent approaches 
its mean field value (/?„ = 0.5) from above p|. Our results for several values 
of the gauge coupling [3 suggest also that the value of 6 increases systemati- 
cally along the critical line with increasing j3, in contrast with the magnetic 
exponent results which are systematically decreasing with j3. 

In spite of the mean field approach for the fermion field, we believe that 
our qualitative picture is realistic. It is in fact hard to imagine that non mean 
field exponents in a mean field approach will become mean field exponents 
after removing the mean field approximation, i.e. that restoring the full 
fiuctuations of the fermion fields would drive back the critical exponents to 
mean field values. 

A numerical analysis of the fermion-gauge-scalar model with compact 
f/(l) gauge symmetry done in 0, has shown the existence of a critical line 
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separating a chirally broken phase from a symmetric phase. In the infinite 
gauge couphng hmit the model is effectively described by the Nambu-Jona 



Lasinio model [ITU] and therefore critical exponents at this point of the criti- 
cal line are gaussian. However strong evidence for non mean field exponents 
has been found in ||^ near the tricritical point separating the second or- 
der line from the first order one. Contrary to our results for the Gauged 
Nambu-Jona Lasinio model, the critical behaviour along the critical line 
of the fermion-gauge-scalar model seems to be well described by the pure 
Nambu-Jona Lasinio model except near the tricritical point. However a slow 
variation of the critical exponents with incresing inverse gauge coupling j3 in 
this model, like the one we have found in the Gauged Nambu-Jona Lasinio 
model, can not be excluded, we believe. 

One important point which deserves further investigation in the (GNJL) 
model is the physical origin of the non mean field behavior. In the (3 ^ oo 
limit of noncompact QED the second contribution to the chiral condensate 
behaves like log m and this result will probably be true also in perturba- 
tion theory. Therefore some important role of gauge configurations topolog- 
ically non equivalent to the free field configuration is suggested by our result 

LU <2. 

We thank CICYT (Spain) - INFN (Italy) collaboration for partial finan- 
cial support to this work. 
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Figure captions 



Figure 1. Phase diagram of the GNJL model in the (3, G plane. 

Figure 2. Logarithm of the nonlinear susceptibility against the logarithm 
of the lattice size for lattice sizes 4,6,8, 10 and (3 = 0.237. 

Figure 3. Inverse generalized nonlinear susceptibility against the inverse 
lattice volume at /3 = 0.237. 

Figure 4. Infinite volume limit of the generalized nonlinear susceptibility 
against q at (3 — 0.237. 
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